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1.0 SUMMARY 


This document identifies the Level C requirements for the Shuttle 
Mission Control Center (SMCC) orbital guidance software. The 
requirements are formulated to provide a functional simulation of 
the onboard Powered Explicit Guidance (PEG) software as required 
by various processors in the SMCC [e.g., Mission Plan Table (MPT), 
Powered Flight Numerical Integrator (PFNI), and Deorbit Planning 
Processor (DPP)]. These Level C guidance software requirements are 
consistent with the Level B requirements of those SMCC processors 
which require a guidance simulation (see References 1, 2, and 3) 
and at a minimum reflect the current status (July 1976) of the 
onboard PEG software as detailed in the preliminary Level C 
Functional Subsystem Software Requirements (FSSR). (See Reference 4.) 

This document presents the formulation of Level C requirements for 
the SMCC guidance software. Section 3.0 presents detailed require- 
ments for a PEG supervisor which controls all input/output interfaces 
with other SMCC processors and determines which PEG mode is to be 
utilized. Sections 4.0 and 5.0 present the detailed description 
(i.e., logic flow diagrams, mathematical equations, etc.) of the 
two guidance modes (7 and 4, respectively) for which these Level C 
requirements have been formulated. Sections 6.0 and 7.0 define 
•functions required for proper execution of the guidance software. 
Section- 6.0 defines the requirements for a navigation function that 
is used in the prediction logic of PEG mode 4. This function is 
extracted from the current navigation FSSR (Reference 5). The 



function PRECISE PREDICT as defined in Section 6.0 may be used 
in the SMCC guidance software for prediction of the burnout state 
until such time that a more concise replacement function is 
defined. Section 7.0 contains routine LTVCON which is used in 
determining a burn/coast-to-target guidance solution. The LTVCON 
routine is defined separately since it will perform a function re- 
quired within the SMCC by users other than PEG. (See Reference 3.) 
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2.0 INTRODUCTION 


The generation of accurate mission profiles (both nominal and abort) 
in the SMCC necessitates a functional simulation of the onboard 
guidance software. As specified in the Level B documents for the 
MPT, PFNI and DPP (References 1, 2, and 3), the requirement exists 
for the simulation of PEG modes 4 and 7 within the SMCC. PEG-4 
(linear terminal velocity constraint guidance) is the primary 
guidance algorithm used during nominal and abort ascent maneuvers 
involving the Orbital Maneuvering System (OMS) engines and for all 
deorbit burns. PEG-7 (guidance for external delta velocity maneuvers) 
is the guidance algorithm for all on-orbit maneuvers and serves as 
the back-up guidance algorithm for all PEG-4 maneuvers. 

The orbital guidance requirements of this memo have been defined by 
McDonnell Douglas Technical Services Company under the Space Shuttle 
Engineering and Operations Support (SSEOS) contract. The requirement 
definitions were developed for the Guidance and Dynamics Branch (GDB) 
of the Mission Planning and Analysis Division (MPAD) in support of 
the- GDB response to Reference 6. 
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3.0 GUIDANCE SUPERVISOR 


This section defines the requirements for a routine (PEGSUP) which 
performs the supervisory function in the orbital guidance software. 
This routine is called by the various SMCC processors and provides 
all necessary input/output and data interfaces needed to provide a 
functional simulation of the onboard guidance software. 


3.1 Detailed Description 

PEGSUP is called using the inputs as detailed in Reference 2. Input 
requirements for PEGSUP are contained in Table I. The table contain 
a list of required inputs for the orbital guidance simulations (e.g. 
target vectors, gravitational constants, propulsion and mass char- 
acteristics). Using the indicated inputs, the function of PEGSUP 
may be summarized by the following sequence; 

a. Test first pass flag (K FpAS$ ) and initialize guidance 

request flags and initialization flag as required. 

b. Test the guidance specification flag (KDq^q) to 
determine if PEG-4 or PEG-7 is to be used. If PEG-4 
guidance is desired, determine if the maneuver is a 


nominal, abort, or deorbit and set appropriate flag, 

c. Compute the change in sensed velocity (AV<.) and update 
the guidance position and velocity vectors (Rq D > Vqq) 


d\- 


and time (Tg D ). 

Determine if active steering computations are desired 
(computation of thrust direction Uy). 
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e. Decrement time-to-go (Tqq) and velocity-to-go (Vqq) if 
active guidance has been terminated. 

f. Update the estimated vehicle mass and acceleration (AIR). 

g. Cycle the guidance computations until a converged solution 
is achieved. 

h. Determine if active guidance computations are to be 

4 

terminated. 

3.2 Logic Flow 

The logic flow of PEGSUP is illustrated in Figure 1. Definitions 
of symbols used in the logic flow diagram are included in Table II. 
Included in Table II is an indicator (TYPE) of the function of each 
variable: 

TYPE DEFINITION 

I Parameter is an input variable to 

this routine 

C Parameter is computed in this 

routine 

0 Parameter is an output from this 

routine to the functional guidance 
routines 

GO Parameter is a guidance output 
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4.0 EXTERNAL DELTA- VELOCITY ROUTINE 

This section defines the routine PEG-7 which performs the computa- 
tions necessary to functionally simulate the onboard external delta- 
velocity guidance algorithms. This guidance mode is primarily used 
to perform powered flight guidance during normal on-orbit maneuvers 
and is designated as the backup guidance algorithm to the linear 

terminal velocity constraint mode (PEG-4) during ascent and deorbit 

* 

maneuvers. 

4.1 Detailed Description 

This guidance routine is called by the guidance supervisor PEGSUP 
to simulate the external delta-velocity algorithm. Using the input 
velocity-to-be-gained vector (V^), estimated vehicle acceleration 
(ATR) and total effective exhaust velocity (V^) of the propulsion 
system, PEG-7 determines the time-to-go until thrust termination 
(Tqq) and the resulting thrust direction (unit vector Uy). On the 
first call to PEG-7, the thrust vector turning rate (X) is set 
equal to zero and is not updated in subsequent passes through PEG-7. 
This insures that the resulting thrust direction (Uy) is equivalent 
to the desired thrust direction (X) for PEG-7 operation.' 

4.2 Logic Flow 

The routine PEG-7 is called from the guidance supervisor PEGSUP 
whenever the guidance mode flag U^jQQjr) is equal to 7. The 
detailed logic flow of this routine is shown in Figure 2. The 
symbols used in the PEG-7 flow diagram are defined in Table III. 
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5.0 LINEAR TERMINAL VELOCITY CONSTRAINT GUIDANCE ROUTINE 


This section defines PEG-4 the routine which provides a functional 
simulation of the linear terminal velocity constraint guidance 
algorithm as currently defined in References 4 and 7. This mode 
performs guidance computations necessary to ensure achievement of 
the desired relationship between horizontal and vertical velocity 
components at a given target position vector. As currently planned 
all OMS ascent maneuvers and the deorbit maneuver will be simulated 
utilizing this mode as the primary guidance algorithm. 

5. 1 Detailed Description 

This routine (in conjunction with its utility routines) determines 

the velocity-to-be-gained {\ Q )> time-to-go until thrust termination 

(Tqq) and the resultant thrust direction vector (Uy) required to 

solve the linear terminal velocity constraint guidance problem. 

Having been called by the guidance supervisor (PEGSUP) with the 

proper input parameters (Table I), this routine then sequences the 

following utility routines to perform the guidance computations: 

Utility Routine Function Performed 

1NIT Parameter initialization (Section 5.1.1) 

UPVGO Update velocity-to-go vector (V rn ) 

(Section 5.1.2) 

TGO Calculate time-to-go (Tqq) until thrust ' 

termination (Section 5.1.3) 

INTEG Calculate thrust integrals (Section 5.1.4) 
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Utility Routine 

Function Performed 

TURNR 

Calculate turning rate vector (A) 
(Section 5.1.5) 

PREDICTOR 

Predict position and velocity vectors 
at cutoff (Section 5.1.6) 

CORRECTOR 

Compute Tqq and test for convergence of- 
of guidance solution (Section 5.1.7) 

PRECISE PREDICT 

Perform a state vector prediction 
function (Section 6.0) 

LTVCON 

Compute the velocity required at an 
initial point in order to intercept a 
target position with a specified linear - 
relationship between the terminal radial 
and horizontal velocity components 
(Section 7.0) 


Upon first entry into PEG-4, INIT is called for parameter initializa- 
tion. On all subsequent passes, UPVGO is called to update and 
the current guidance phase number. If V^q is zero, PEG-4 is exited 

J 

(no guidance computations are required). After determining Tqq, the 

\ T 

desired thrust vector (X), and the turning rate vector (A), a test 
is made to check for unit thrust direction (Uj) calculation. If 
active steering is enabled ^STEER~^ anc * ^ as com Puted a 

converged solution with no significant differences between the current 
and previous predicted burnout times, Uy is computed. Otherwise this 
calculation is bypassed. Following thrust direction calculations, 
burnout position and velocity-to-go are determined. 

The linear terminal velocity constraint logic (PEG-4) is called from 
PEGSUP to perform mode 4 guidance computations. The detailed logic 
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flow of PEG-4 is shown in Figure 3. Definitions of symbols used in 
PEG-4 and its utility routines are found in Table IV. 

5.1.1 Utility Routine INIT 

This routine is called only on the first pass through the mode 4 
guidance computations. INIT sets the predicted burnout state 
(Rp, Vp) equal to the current' state and calls utility routine 
CORRECTOR to determine an initial V^. INIT also initializes the 
current guidance phase indicator to one. 

A detailed logic flow of INIT is shown in Figure 4. Definitions of 
symbols used in this routine are found in Table IV. 

5.1.2 Utility Routine UPVGO 

On all passes through PEG-4 subsequent to the first pass, utility 
routine UPVGO is called to update the velocity-to-be-gained before 
thrust termination (Vgg). V^q is decremented by the computed change 
in accumulated sensed velocity (AV^) as determined in PEGSUP. UPVGO 
also tests to determine when the current guidance phase is terminated 
and increments the phase indicator ( ^phASe) 9 in '* tial ' 5 * zes the burn 
time of the new phase (Tg(K pHASE )) and initializes guidance phase . 
time-to-go ( t qo/\ ^ K PHASE^ * A detailed lo 9^ c f1ow of routine 

UPVGO is shown in Figure 5. Definitions of symbols used in the 
flow diagram are found in Table IV. 

5.1.3 Utility Routine TGO 

Utility routine TGO is called by PEG-4 to compute time-to-go until 
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thrust termination (Tgg) and to update predicted burnout time (T ). 

A test is made to determine if multi-phase logic (N>1 ) is to be 
utilized to compute Tgg. If single phase logic is specified, Tgg 
is computed as a function of the vehicle estimated acceleration 
(AIR), total exhaust velocity (V^) and magnitude of the velocity- 
to-go (V G0MAG ). If multi -phase logic is utilized, estimated burn 
times of all phases are computed and summed to determine the total 
T g0 for the maneuver. Following the Tgg computation, predicted burn- 
out time (T ) is computed by summing current time and Tgg. It should 
be noted that the multi -phase software requirements specified in 
Figure 6 are limited by: a) multi-phase implies two phases only 

b) multi -phase is designed for an OMS/RCS 

parallel burn followed by an OMS burn. 

« 

Detailed logic flow of routine TGO is shown in Figure 6. Definitions 
of symbols used in the logic flow diagrams are found in Table IV. 

5.1.4 Utility Routine INTEG 

INTEG is a utility routine called by PEG-4 to compute a set of inputs 
for determining the first and second order thrust integrals. The . 
integral calculations assume a constant force magnitude and are time 
integrals of force over mass. The integrals are used in the thrust 
vector turning rate calculations and the prediction equations ot 
guidance. The integrals are computed for each phase and summed over 
the total number of guidance phases. 

As currently defined in Reference 7, the deorbit and AOA OMS-2 
guidance algorithms consist of single phase logic, utilizing only 
.first order thrust integrals. 
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The orbital guidance requirements specified in this report are 
consistent with the reference. The computation of second order 
integrals and multi-phase logic are deleted depending on a test 
which is incorporated into INTEG to determine if deorbit/AOA OMS-2 
guidance is desired. 


After determining the thrust integrals, INTEG computes a reference 
time (T^) used in the calculation of the resultant thrust direction 

(u T ). 

Detailed logic flow of INTEG is illustrated in Figure 7. Definitions 
of symbols used in the logic flow diagram are given in Table IV. 


5.1.5 Utility Routine TURNR 

Routine TURNR is called by PEG-4 to determine .thrust vectors X and 
X. . The Jesired thrust vector (X) is computed as a unit vector in 


the direction of V^q. If the simulated maneuver is a deorbit or 
AOA OMS-2, the turning rate vector (X) is computed as a function of 
the primer rate (X^), thrust vector (X) and a unit vector (Uy) 
which is normal to the desired orbital plane. If the maneuver is ‘ 
not a deorbit or AOA OMS-2, X is computed as a function of the 
range-to-go vector (Rqq), X and the thrust integrals. A test is 
made to insure that X does not result in a violation of small angle 
approximations used in development of' the guidance equations. 


Detailed logic flow of TURNR is illustrated in Figure 8. Definitions 
of symbols used in the logic flow diagram are found in Table IV. 
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5.1.6 Utility Routine PREDICTOR 

Routine "PREDICTOR is called by PE6-4 to predict the burnout state 
vector. To perform this function, the first and second integrals 
of thrust acceleration and gravity must be determined. The function 
of PREDICTOR and its interface with a navigation-defined state 
propagation routine are described in this section. 

5. 1.6.1 Determination of Thrust Integrals 

The thrust time integrals (VjhrusT 5 Rjhrust^ are com P ute ^ as 
functions of the reference thrust vectors (A., A) and the force over 
mass integrals of INTEG for all simulated maneuvers except deorbit 
and AOA OMS-2. As specified in Section 5.1.4, Reference 7 indicates 
that the second order integrals of force over mass are negated for 
these maneuvers. Hence, the thrust time integrals are determined 
by setting \T THR y ST equal to V^q and ec l ua ^ 1:0 ^gq ^ or c5eorf3 '’ t 

and AOA OMS-2. 

5 . 1 . 6 . 2 Determination of Gravity .Acceleration Integrals 
Integrals of gravity acceleration are approximated by propagation of 
the state vector over a coasting arc which remains close to the 
powered trajectory. The initial state vector (R^ 5 V C1^ input 

to the state propagation routine are computed as a function of the 
current -guidance state vector (Rqq* v gd^ 5 thrus1: integrals 
^THRUST 5 ^THRUST^’ anc ^ ^GO* Curren t onboard guidance software 
specifications (Reference 4) indicate that the state propagation is 
to be performed by a navigation function. The logic for this 
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navigation function has been extracted from the navigation FSSR 
(Reference 5) and is defined in this memo as routine PRECISE 
PREDICT. A description of PRECISE PREDICT is given in Section 6.0. 

As part of the inputs to PRECISE PREDICT, PREDICTOR determines a 
maximum value for integration step size (AT^g). Step size AT^yg 
is set equal to the middle value of ATw*w, T G0 /N SEG and AT MIN' 
in PREDICTOR. 

5. 1.6. 3 PREDICTOR/PRECISE PREDICT Interface 
Targets for a PEG-4 maneuver are based on a burn-coast sequence 
(input target values are for the end coast phase). The burn arc 
of the PEG-4 solution includes oblateness effects (0^ f 0) in the 
gravity potential model. The coast arc within PEG-4 is calculated 
using a spherical earth model. Thus it is necessary to bias 
guidance input targets in order to achieve the desired terminal 
state vector. The desired terminal conditions (Ry, Cp C 2 ) are 
biased (Ry + ARy, C-j + ACy ) to account for oblateness effects (J 2 ) 
which will be experienced during the actual coast. The target 
biases (ARy and AC-j) are inputs to the guidance software as defined 
in this paper. 

The target biases are computed based on a given burn arc with no 
engine failure; nominal burn arc is assumed to be accurately pre- 
dicted.* - However, to account for off-nominal situations (predicted 
burn arc does not equal given burn arc), the position at cutoff 
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(Rp, Vp) must be adjusted for the oblateness effects. This neces- 
sitates two calls to PRECISE PREDICT. The first call propagates 


Rq-j- and Vqi (with effect included) over a time interval 
equivalent to T^q minus the difference between the targeted coast 
time. (AT C q£ ST and the actual coast time (^TqoAST^* 

second PRECISE PREDICT call propagates the state resulting from 
the first call (R- INT > Vj^j) over an interval equivalent to' ' 

AT C0AST (REF) minus AT C0AST (-Without J 2 effect). The final state 
vector is defined as Rq 2 ’ ^C2' 


Having accomplished the state extrapolation, the gravity biases 
(^GRAV’ ¥ GRAV ^ are com P u1:ecJ and are used to adjust the predicted 
burnout state (R , V ). 

r r 


5.1.6. 4 PREDICTOR Logic Flow 

Detailed logic flow of routine PREDICTOR is illustrated in Figure 9 
Definitions of symbols used in the logic flow diagram are found 
in Table IV. 


5.1.7 Utility Routine CORRECTOR 

CORRECTOR is called by PEG-4 to set the desired position vector 
(R D ), to determine the desired velocity (Y Q ) at burnout, to correct 
VgQ for any desired fuel depletion burns and to perform the check 
on guidance convergence. Desired velocity (V" D ) is determined by 
routine* LTV CON. LTVCON computes V Q such that the target position 
(at end coast) is intercepted with a desired relationship between 
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horizontal and vertical velocity. Before the LTVCON call is executed, 
corrections are made to the target position (Ry). If the simulated 
maneuver is an abort or deorbit, Ry is constrained to be in the 
predicted cutoff plane. If the maneuver is an AOA OMS-2 or deorbit 
maneuver, the magnitude of Ry is adjusted to be consistent with the 
desired geodetic altitude (Hygy) above the Fischer ellipsoid. 


After determining Vp, a correction is applied to V GQ using a scalar 
damping factor (RHO^g) and miss velocity (difference between 7 and 
Vq). If fuel depletion burns are to be .performed, Vqq is. corrected 
to include the out-of-plane velocity accumulated by the depletion 
. burn. After correcting Vgg, a convergence check is made to determine 
if miss .velocity (^ 55 ) is within acceptable limits. 


Detailed logic flow of routine CORRECTOR is illustrated by Figure 10. 
Definitions of symbols used in the logic flow diagram are given in 
Table IV. 
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6.0 STATE PREDICTION ROUTINE 

This section defines routine PRECISE PREDICT. The routine is called 
by PREDICTOR to perform the state extrapolation function. The 
extrapolation is required to determine gravitational effects to 
be used- in the prediction of burnout position and velocity vectors 
(Rp, Vp). PRECISE PREDICT is developed from a modification of the 
current "Super G" algorithm. (See Reference 5.) 

6.1 Detailed Description 

The function of PRECISE PREDICT is to calculate an extrapolation of 
a state vector. To accomplish this the following inputs are 
required: 

a) initial state vector - (Rp Vp. Tj) 

b) ' time for computation of final state (Tp) 

c) upper bound on integration step size (AT^) 

d) indicators for gravity, drag and- vent models (GMD, GMO, DM, VM) 

Using the above inputs, PRECISE PREDICT determines an integration 
interval for use in the calculations and. the initial acceleration. 
The routine then: 

a) advances current time (T) 

b) integrates current radius (Rp) and velocity (Vp) 

c) updates acceleration 

d) corrects radius by the "Super G" algorithm. 
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This process is repeated iteratively until time Tp. The last pass 
values of and V r are the state vector at the final time. 

Sections 6.1.1 and 6.1.2 contain descriptions of two utility 
routines (ACCEL and ACCEL GRAV) which are. used by PRECISE PREDICT. 

Detailed logic flow of routine PRECISE PREDICT is illustrated by 
Figure 11. Definitions of symbols used in PRECISE PREDICT-and its 
utility routines are contained in Table V. 

6.1.1 Utility Routine ACCEL 

ACCEL is a function routine. It is called from PRECISE PREDICT 
to compute gravitational and non-gravitational (drag, venting) 
accelerations. The SMCC requirements for drag and venting models 
(within the guidance link) are unclear at this- time. (See Reference 
7.) These models are controlled by input flags DM and VM and are 
not called if DM = 0, VM = 0. Because the requirement for these 
models is not clearly defined at this time, only the logic for 
accessing drag and venting models is included in this report. 

ACCEL is the supervisor for the acceleration calculations. It computes 
acceleration due to a central body force, controls the calling of a 
higher order gravitational model (ACCEL GRAV), controls calling of 
drag and vent models (ACCEL DRAG and VENT) , and computes the total 
acceleration. In addition, it issues calls for vehicle attitude (ATT) 
and atmospheric density (DENSITY). Vehicle attitude is required by 
VENT and ACCEL DRAG. An attitude control flag (ATM), is set in the 
prediction routine. The flag is consistent with the Reference 5 logic 
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in that the referenced attitude determination logic requires an 
indicator for propagation (ATM = 0) or prediction (ATM =1). 

Logic flow for routines ACCEL and ACCEL GRAV is provided in this 
report.. The vehicle attitude (ATT), vent (VENT), density (DENSITY) 
and drag (ACCEL DRAG) models are not included. The outputs re- 
quired from these routines are specified in Figure 12 and Table V. 

Detailed logic flow of the function ACCEL is illustrated by Figure 
12. Definitions of symbols used in the logic flow diagram are 
listed in Table V. 

6.1.2 Utility 'Routine ACCEL GRAV 

This section contains a brief description of utility routine ACCEL 
GRAV. The logic flow of ACCEL GRAV is illustrated in Figure 13. 

The logic for the model is taken from Reference 5 and is based on 
Reference 8. This model determines gravitational effects resulting 
from the oblateness of the Earth. The model is controlled by input , 
flag GMD (GMD f 0: compute oblateness effects). Input flag GMO 

indicates the order of calculation for the model. The inputs are 
restricted in that GMO _< GMD with a maximum value of 8 for GMD. The 
limits which are imposed on the arrays of ACCEL GRAV are specified 
in Table V. - 

The logic contained in Figure 13 for determining gravitational accelera- 
tion includes both sectorial and zonal effects. The acceleration 
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calculations are performed in earth-fixed coordinates and rotated into 
the mean of 1950 (M50) system for use by ACCEL. A detailed descrip- 
tion of the logic of ACCEL GRAV is contained in Reference 8. Table V 
contains definitions of the symbols used in the Figure 13 logic flow. 
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7.0 LINEAR TERMINAL VELOCITY CONSTRAINT ROUTINE 


This section describes routine LTVCON which is used to determine 
the velocity required at an initial point ;Rg) to intercept a 
target position (fU) such that a specified linear relationship be- 
tween the terminal radial and horizontal velocity components is 
satisfied. This linear relationship is specified by: 


v radial c i + c 2 * Horizontal 

The solution involves a set of two-body differential equations of 
motion for a spacecraft center of mass in an inverse square central 
force field. 


7.1 Detailed Description 

Within the logic flow as specified for calculating required 
velocity, LTVCON is an executive routine. LTVCON controls 
the calculation of velocity-required through the calling of 
three utility routines. The routines and their function are as 


specified: 
UTILITY ROUTINE 
INITIAL 

REQVEL 


TRANSTIME 


FUNCTION ' 

Compute unit vector specifying transfer plane. 
Determine transfer triangle and constants for 
normalizing velocity and transfer time. 

Compute velocity vectors at Rg and R, such that 

terminal velocity satisfies the desired linear 
constraint. 

Compute time required to traverse the transfer 
angle given a specified initial velocity (Vg). 
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The logic flow of LTVCON is illustrated by Figure 14. Table VI 
contains definitions of the variables used in the entire LTVCON 
link including the specified utility routines. 

7.1.1 Utility Routine INITIAL 

Utility routine INITIAL performs an initialization operation. It 
computes the magnitudes and unit vectors for input position vectors 
R 0 and R-p The routine also computes sines (S6) and cosines (C9) 
of the transfer angle and sets the constants (Rqjrq arid V^q) 
required to normalize transfer velocity and time. 

Figure 15 contains an illustration of the logic flow of INITIAL. 
Definitions of the symbols used in the figure are 'contained in 
Table VI. 

7-1-2 Utility Routine REQVEL 

REQVEL computes the velocity vector (Vg) required at an initial 
position (R ) in order to achieve the required velocity vector (V-j) * 

at a final position (R^). Since the terminal radial velocity component 
(V Rn ) is specified as a linear function of the terminal horizontal 
velocity component (V^), a quadratic equation in V H1 exists as 
follows: 

A * V m 2 + B * V R1 + C = 0, 

where the coefficients A, B, and C are functions of the- position - 

magnitudes • (R n , R. ), the transfer angle, and the constraint 
U MAG 'MAG 
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coefficients (C-j , C 2 ). The radical of the equation (D) is computed 
as: 

' D = B 2 - 4 (A * C) 

If D is less than zero, no solution to the problem exists and an 
error flag is set and a message displayed. At this point LTVCON and 
the calling routine should be exited. If D is greater than or equal 
to zero, the velocities V Q and V-j are computed such that the desired 
velocity constraint condition is fulfilled. 

The logic flow of REQVEL is illustrated in Figure 16. Definitions 
of the symbols used in the figure are contained in Table VI. 

7.1.3 Utility Routine TRANSTIME 

Routine TRANSTIME computes the time required to traverse from 
Wq to R”-j. To determine the free-fall time required to traverse the 
transfer angle, an indicator (U) is computed to determine if a 
solution exists. If U is less than or equal to a -1 , the REQVEL 
solution results in a hyperbolic trajectory. In this case an error 
flag is set, an error message displayed and LTVCON and its calling 
routine should be exited. If U is greater than -1, transfer time 
(AT) is computed. The computation of AT requires the use of a- 
function QM. This function routine computes 

QM = W + (1 - WW) (5 - (F) ( WW ) ) / 1 5 

where F is a hypergeometric function which is the transcendental 
portion of the normalized transfer time. 
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Figure 17 contains an illustration of the logic flow of TRANSTIME. 
Figure 18 contains the logic flow of function QM. The symbols 
used ir. the figures are defined in Table VI. 
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8.0 CONCLUSIONS 


The logic and logic flow given in this report are developed to 
satisfy the Level C orbital guidance software requirements for 
the SMCC. The resulting functions are to provide guidance simu- 
lation capability within the SMCC. The logic and flow presented, 
in this report are modeled after the current guidance FSSR but are 
not restricted to being identical. The routines used to perform 
guidance functions are as depicted in Figure 19. 

Some of the functions defined in this report may be performed by 
other disciplines within the SMCC. Orbital guidance is only a 
"user" of these functions (e.g. the gravity model). Some of these 
functions are defined in this report for completeness. In some 
instances the capability (or requirement) for the functions (e.g. 
vent model and drag model) and its interface with the models of 
this report are defined. . r , _ 

As of this date it is understood that here is no requirement for 
simulating the drag and vent models within the orbital guidance 
function of the SMCC. The logic for these two routines should be 
removed from these requirements. 
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TABLE I 

- DEFINITIONS' OF 

INPUTS TO 

SMCC GUIDANCE SOFTWARE 

SYMBOL 

UNITS 

DIMENSION 

USER 

ROUTINES* 

DEFINITION 

C 1 

FT/SEC 

1 

C 

Target velocity constraint coefficient in the 
equation V RmAl = C ] + C £ * V HQRIZ 

C 2 

NA 

1 

C 

Target velocity constraint coefficient ‘in the 
equation V RADIAl = C-, + * V HQRIZ 

C L 

NA 

L-l, 35 

' C 

Constants which account for tesseral harmonics 
in gravity potential model 

[EF TO M50(T)] 

NA 

3,3 

C 

Transformation- matrix from earth-fixed to mean 
of 1950 (M50) coordinates 

ELLIPT 

NA 

1 

C 

Earth ell ipti city constant 

h tgt 

FT 

1 

c 

Geodetic altitude of target above Fischer el 1 ipsoid 

f oms 

LB 

1 

A/C 

Nominal OMS vacuum thrust level 

f rcs 

LB 

1 

A/C 

Nominal RCS vacuum thrust level 

^GUID 

NA 

1 

A . 

Guidance control option 
= 1; nominal PEG-4 maneuver 
= 2; PEG-7 maneuver 
= 3, 4, 5; manual steering modes 
= 6; deorbit maneuver 
= 7; abort maneuver (AOA or ATO) 
= 8; AOA OMS-2 maneuver 

k fpass 

NA 

1 

A 

First pass flag 
= 0; subsequent passes 


= 1; first pass through guidance 


*A = PEGSUP; 3 = PEG-7; C = PEG-4 




TABLE 

I - DEFINITIONS 

OF INPUTS TO 
(CONTINUED) 

SMCC GUIDANCE SOFTWARI 

SYMBOL 

UNITS 

DIMENSION 

t 

USER 

, ROUTINES* 

DEFINITION 

k fueld 

NA 

1 

C 

Fuel depletion flag 
= 0; no fuel wasting 
= -1; waste right 
= +1; waste left 

Sm 

NA 

1 

c 

Drag model flag used in PRECISE PREDICT 
= 0; no drag 

= 1; compute drag acceleration 

■Stm 

NA 

1 

c 

Vent model flag used in PRECISE PREDICT 
= 0; no venting 

= 1; compute vent accelerations 

MASS 

SLUGS 

1 

A/C 

Mass of orbiter at start of maneuver. 

M rr • 
tta OMS 

SLUGS 

1 

A/C 

Mass of OMS propellant burned during a fuel 
depletion maneuver to achieve the desired 
orbiter mass 

^OMS 

SLUGS/SEC 

1 

A/C 

Nominal OMS mass flow rate 

\cs 

SLUGS/SEC 

1 

A/C 

Nominal RCS mass flow rate 

N MAX ' • 

NA 

1 1 

A 

Maximum number of guidance iteration cycles 
allowed ■ ■ 

n phase 

NA 

1 

A 

Number of guidance thrust phases 

n seg 

SEC 

1 

C 

Number of integration time steps used for 
gravity prediction 

R 

FT 

3 

A 

Current vehicle inertial position vector (M50) 

r eq 

FT 

1 

C 

Earth equatorial radius 

*A = PEGSUP; 

B = PEG-7* C = PEG-4 






TABLE I - DEFINITIONS OF INPUTS TO SMCC GUIDANCE SOFTWARE 

(CONTINUED) 


USER 


SYMBOL 

UNITS 

DIMENSION 

ROUTINES* 

DEFINITION 

r egrav 

FT 

* 

1 

C 

Radius of earth used 

for gravity computations 
« 

r t 

FT 

3 

C 

target radius vector (M50) 

S, 

L 

NA 

L=l, 35 

c. 

Constants which account for sectorial harmonics 
in gravity potential model 

S MISS 

NA 

1 

C 

Decimal fraction of 
criterion for I^vjiss 

V G0 | used as a convergence 

T 

SEC 

1 

A 

Current time (GMT) 


T 

SEC 

1 

A 

Time before thrust termination at which active 

go min 


guidance computations are inhibited 

r iG 

SEC 

1 ' 

A 

Maneuver ignition time 

T RCSOFF 

SEC 

1 

C 

RCS cutoff time during AOA OMS-1 parallel OMS/ 
RCS maneuver 

U’gp 

NA 

3 

c 

Unit vector in direction of the Earth's axis 
of rotation 

V 

FT/SEC 

3 

A • 

Current vehicle velocity vector (M50) 

V s 

FT/SEC 

3 

A 

Current vehicle sensed velocity vector 

V GO 

FT/SEC 

3 

B 

Required velocity-to 

•rbe-gained vector (PEG-7 ) 

XLD XZ 

NA 

1 

C 

Fractional multiplier of circular orbital rate 
used in computation of primer rate (X^) 

ZONAL j 

*A = PEGSUP; B 

NA 

= PEG-7-' C = PEG-4 

1=1 ,'8 

C 

Constants which account for zonal harmonics in 
gravity acceleration model 



TABLE I - DEFINITIONS OF INPUTS TO SMCC GUIDANCE SOFTWARE 

(CONTINUED) 


USER 


SYMBOL 

UNITS 

DIMENSION 

ROUTINES* 

DEFINITION 

AC-i 

FT/SEC 

1 

C 

Offset applied to true target velocity con- 



straint C-| to obtain a conic constraint 

ARj 

FT 

3 

A 

. 

■ Offset applied to true target position vector 
Rj to achieve a conic target 

AT COAST(REF) 

SEC 

1 

c 

Reference coast time used in prediction logic 

at cutoff 

SEC 

1 

A/C 

Effective thrust tail -off time 

at max 

SEC 

1 ' 

. C 

Maximum integration time step for gravity 
prediction 

aT min 

SEC 

V 

C 

Minimum integration time step for gravity, 
prediction 

e STEER 

SEC • 

1 

C 

Convergence criterion applied to predicted _ 
burnout time used to inhibit computation of \ij 

^MAX 

RAD 

1 

C 

Small angle approximation constraint 

^ EARTH 

ft 3 /sec 2 

1 

C 

Earth gravitational constant. 


*A = PEGSUP; B = PEG-7; C = PEG-4 



FLOW CHART 
SYMBOL 

ATR 

f oms 

f rcs 

FT 

I FLAG 


K ABORT 

K A0MS2 

!< conv • • 
Cutoff 

^OEORB 

*C = COMPUTED; I 


TABLE II - DEFINITIONS OF SYMBOLS USED IN PEGSUP 


UNITS 

. DIMENSION 

TYPE* 

FT/SEC 2 

1 

C/0 

LB 

1 

I 

LB 

1 

I • 

LB 

1 

C 

NA 

1 

C/0 


NA 

1 

C/O 

NA 

1 

C/0 

NA 

1 * 

c 

NA ' 

1 

•c 


DEFINITION 

Current estimated thrust acceleration 

Nominal OMS vacuum thrust level 

Nominal RCS vacuum thrust level 

Total vehicle vacuum thrust 

PEG error flag 
= 0; no error 

= 1 ; no physical solution (from LTVCON) 

= 2; hyperbolic trajectory (from LTVCON) 

= 3; Vg g equals zero 

Abort guidance request flag 
= 0; no abort simulation 
= 1; abort simulation (AOA or ATO) 

OMS-2 burn guidance flag (used only when 
K AB0RT = 1 ^ 

= 0; AOA OMS-2 guidance not desired 
= 1; perform AOA OMS-2 guidance 

Convergence flag 

= 0; PEG solution has not converged 
= 1; PEG solution has converged 

Active guidance termination flag 
= 0; do not perform guidance cutoff logic 
= 1; perform guidance cutoff logic 

Deorbit guidance request flag . 

= 0; not deorbit maneuver 
- 1; deorbit maneuver 


NA 1 C/0 

INPUT; 0 = OUTPUT TO PEG-4 OR PEG-7; GO = GUIDANCE OUTPUTS 



TABLE II - DEFINITIONS OF SYMBOLS USED IN PEGSUP 

(CONTINUED) 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

TYPE* 

DEFINITION 

K FPASS 

. NA 

1 

I 

First pass flag 

= 0; subsequent passes 

= 1; first pass through guidance 

k init 

.NA 

1 

C/O 

Initial pass flag (set to 1 on initial 
pass) 

’Siode 

NA 

1 

C 

Guidance mode option 
« 4; use PEG-4 
=7; use PEG-7 

k steer 

NA 

1 

C/O 

Active steering flag_ 
= 0; do not compute U~ 
- 1 ; compute 

■^GUID 

NA 

1 

I 

Guidance control option 
« 1; nominal OMS burn guidance (PEG-4) 

= 2; external A V guidance (PEG-7) 

= 3 3 4 s 5; manual steering modes 
= 6; deorbit maneuver (PEG-4) 

= 7; abort maneuever (AOA or ATO) (PEG-4) 
= 8; AOA OMS-2 maneuver’ (PEG-4) 

mss 

SLUGS 

1 

, I/C 

Estimated mass of orbiter 

f>1 BO 

SLUGS 

1 

c 

Desired mass of orbiter at termination of 
fuel depletion maneuver 

M rr 

Ub OMS 

SLUGS 

1 

I 

Mass of OMS propellant burned during a 
fuel depletion maneuver to achieve the 
desired orbiter mass 

*C = COMPUTED; 

I ** INPUT; 0 = 

OUTPUT TO PEG-4 OR 

PEG-7; GO = 

GUIDANCE OUTPUTS 





TABLE II - DEFINITIONS OF SYMBOLS 

USED IN PEGSUP 


• 


(CONTINUED) 


FLOW CHART 




DEFINITION 

SYMBOL 

UNITS 

DIMENSION 

•TYPE* 

%S 

SLUGS/ SEC 

1 

I 

Nominal OMS mass flow rate 

^RCS 

SLUGS/SEC 

1 

I 

Nominal RCS mass flow rate 

M 

SLUGS/SEC 

1 

C 

Total vehicle mass flow rate 

N„.,, 

NA 

1 

C 

Number of calls to PEG-4 or PEG-7 per 

CYC 

_ - - — - 


guidance cycle 

N 

NA 

1 

I 

Maximum number of guidance iteration cycles 

MAX 



si lowed 

to w 

M "phase 

NA 

1 

I 

Number of guidance cycles 

R 

FT 

3 

I 

Current vehicle position vector input to 
guidance 

^GD 

Ft . 

3 

C/O 

Vehicle position vector used by guidance 

T 

SEC 

I 

I • 

Current time 

t cutoff 

SEC 

1 

c 

Current time of engine cutoff 

T ■ 

• *GD 

SEC 

I 

c 

Current guidance time 

‘'GDP 

SEC 

1 

c 

Previous value of guidance time 

t go 

SEC ’ 

1 

C/GO 

Time-to-go to thrust cutoff 

Tgo min 

SEC 

1 

I 

Value of T gq when active guidance is inhibited 

t ig 

SEC 

1 

I 

Maneuver ignition time 

*C = COMPUTED; 

I = INPUT; 0 = OUTPUT TO PEG-4 OR 

PEG-7; GO = GUIDANCE OUTPUTS 



TABLE II - DEFINITIONS OF SYMBOLS USED IN PEGSUP 

(CONTINUED) 

>• \ 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

TYPE* 

DEFINITION 

u T 

NA 

3 

GO 

Unit vector in the thrust direction 

l . 

¥ 

FT/SEC 

3 

I 

Current vehicle velocity vector input to 
guidance 

V EX 

FT/SEC 

1 

C/O 

Total effective exhaust velocity 

^GD 

FT/ SEC 

3 

C/O 

Vehicle velocity vector used by guidance 

^GO 

FT/SEC 

3 

I/C/O/GO 

Velocity-to-be gained vector 

T s 

FT/SEC 

3 

I 

Current sensed velocity 

h? 

FT/ SEC 

3 

c 

Previous pass sensed velocity 

at cutoff 

• SEC 

1 

I 

Effective thrust tail -off time 

at gd 

SEC 

1 

c . 

Difference between current and previous 
values of T^ 

AV $ ' 

FT/SEC 

3 

• C/O 

Change in accumulated sensed velocity from 
previous value 

AV<. 

b MAG 

FT/SEC 

1 ' 

C 

Magnitude of AV<. vector 
*- 

*C = COMPUTED.; I = 

INPUT; 0 = 1 

OUTPUT TO PEG-4 OR PEG-7 

GO = GUIDANCE OUTPUTS 



TABLE III - DEFINITIONS OF SYMBOLS USED IN PEG-7 


FLOW CHART 
• SYMBOL 

♦ 

UNITS 

DIMENSION 

TYPE* 

DEFINITION 

ATR 

FT/SEC 2 

1 

I 

Current estimated acceleration 

k conv 

: • na 

1 

C 

Convergence flag 

= 0; PEG solution has not converged 
= 1; PEG solution has converged 
KcONV ~ for PEG_7 guidance 

k init 

NA 

1 

I 

Initialization request flag 
» 0; no initialization required 
= 1; perform initialization 

• K 3TEER 

NA 

1 

I 

Active steering flag_ 
= 0; do not compute Uy 
= 1 ; compute Uy 

t go 

SEC 

1 

■ C/0 

Time-to-go to thrust cutoff 


SEC 

1 

c 

Reference time for thrust vectors X and X 

V • 

NA 

3 ’ 

C/0 

Unit vector in the resultant thrust direction 

V E X 

FT/SEC 

1 

I 

Total effective exhaust velocity 

? G0 

FT/ SEC 

3 

I/C/O 

Velocity-to-be-gained vector 

V GOMAG 

FT/ SEC 

1 

C 

Magnitude of V^q 

AV S 

FT/SEC 

3 ' 

I 

Change in accumulated sensed velocity 

X 

NA 

3 

c 

Unit vector in direction of desired thrust 

X 

NA 

3. 

c 

Thrust vector turning rate 

X 

SEC 

1 

c 

Ratio of mass to mass flow rate 

*C = COMPUTED; 

I = INPUT; 0 = OUTPUT 






TABLE IV - DEFINITIONS OF SYMBOLS USED IN PEG-4 


FLOW CHART 


SYMBOL 


UNITS 

DIMENSION 

TYPE* 

DEFINITION 

ATR 


FT/SEC 2 

1 

I 

Current estimated thrust acceleration 

C 1 

• 

FT/SEC 

1 

I 

Target velocity constraint coefficient in- the equation 
V RADIAL = C 1 + C 2 V HORIZ 

C 2 


NA 

1 

.1 

Target velocity constraint coefficient in the equation 
V RADIAL = C 1 + C 2 V H0RIZ 

ELLIPT 


NA 

1 

I 

Earth ell ipti city constant 

F OMS 


LB 

1 

I 

Nominal OMS vacuum thrust 

H 


FT-SEC 

1 

c 

2 

Time integral from zero to T^q of (FT/M)T 

HA( ) , 


FT-SEC 

n phase 

c 

Segment of H contributed by each guidance 
phase 

h tgt 


FT 

1 

I 

Geodetic altitude of target above Fischer 
ellipsoid 

I FLAG 


NA 

1 

c 

PEG error flag 
= 0; no error 
= 1 ; no physical solution 
= 2; hyperbolic trajectory 

= 3; equals zero 

0 


FT , 

1 

c 

Time integral from zefo to T^q of (FT/M)T 

OA( ) 


FT 

n phase 

c 

Segment of 0 contributed by each guidance phase 

COL 


' SEC 

1 

c 

Burn centroid time 


*C = COMPUTED; I = INPUT; 0 = OUTPUT 



TABLE IV - DEFINITIONS OF SYMBOLS USED IN PEG-4 

(CONTINUED) 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

TYPE* 

DEFINITION 

k abort : 

NA 

1 

I 

Abort guidance request flag 
= 0; no abort simulation 
= 1; abort simulation (AOA or ATO) 

K A0MS2 

NA 

1 

I 

OMS-2 burn guidance request flag (used only 

when k abort = n ' 

= 0; AOA OMS-2 guidance not desired 
= 1; perform AOA OMS-2 guidance 

k conv 

NA 

i 

c 

Convergence flag 

= 0; PEG solution has not converged 
=1; PEG solution has converged 

k deorb 

NA 

* 


I 

Deorbit guidance request flag 
= 0; not deorbit maneuver 
= 1; deorbit maneuver 

k dm 

NA 

i 

I ■ 

Drag model flag for precise predictor 
= 0; do not compute acceleration due to drag 
= 1; compute acceleration due to drag 

k fueld , 

NA 

i 

I 

Fuel depletion flag 
= 0; no fuel wasting 
= +1 ; waste left 
= -1 ; waste right 


*C - COMPUTED; I - INPUT; 0 - OUTPUT 



TABLE IV - DEFINITIONS OF SYMBOLS USED IN PEG-4 

(CONTINUED) 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

TYPE* ' 

DEFINITION 

K INIT 

NA 

1 

I 

Initialization request flag 
= 0; no initialization required • 
= 1; perform initialization 

K J2 

NA 

1 

C 

J2 compensation flag 

= 0; do not include effects in precise 

predictor 

= 1; include ^ effects in precise predictor 

k phase 

NA 

1 

C 

Current guidance thrust phase indicator 

k steer 

NA 

1 

I 

Flag to enable active. steering 
= 0; do not compute Uj 

= 1 ; compute Uy 

k vm 

NA 

1 

I- 

Vent model flag for precise predictor 
= 0; do not compute acceleration due to venting 
= 1; compute acceleration due to venting 

L 

FT/SEC 

1 

c 

Time integral from zero to Tqq of FT/M 

LA( ) 

FT/SEC 

N PHASE • 

c' • 

Segment of L contributed by each guidance 
phase 

MASS 

SLUGS 

1 ' 

I/C • 

Current estimated mass of orbiter 

M 

SLUGS/SEC 

1 

c 

Current mass flow rate 

m bo 

SLUGS 

1 

c 

Desired mass of orbiter at thrust cutoff 


*C = COMPUTED; I = INPUT; 0 = OUTPUT 



TABLE IV - DEFINITIONS OF SYMBOLS USED IN PEG-4 

(CONTINUED) 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

^oms : 

SLUGS/SEC 

1 

N 

.NA 

1 

n seg 

NA 

1 

p 

FT-SEC 2 

1 

PA( ) 

FT-SEC 2 

n phase 

Q 

FT-SEC 

1 

Q / 

FT-SEC 

1 

QA( ) 

FT-SEC 

n phase 

r eq 

FT 

1 

r gox 

FT . 

T 

1 

r goy 

FT 

1 

r goz 

FT 

1 

*C v COMPUTED ;< I = 

INPUT; 0 = OUTPUT 



TYPE* DEFINITION 

I Nominal OMS mass flow rate 


.1 

I 

C 

C 

c 

c 

c 

I 


Number of active guidance thrust phases used 
by guidance (N » N pHASE ) 

Number of integration time steps for gravity 
prediction 

Double time integral from zero to T then zero 
to T qo of (FT/M)T^ 

Segment of P contributed by each guidance 
phase 

Double time integral from zero to T then zero 
to T gq of ( FT/M )T 

Intermediate variable in thrust integral 
computations 

Segment of Q contributed by each guidance 
phase 

Earth equatorial radius 


C Component of Rqq along U x 
C Component of Rqq along Uy 
C Component of R^q along U^ 



00 
' OD 




TABLE IV - DEFINITION OF SYMBOLS 

(CONTINUED) 

USED IN PEG-4 

FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

TYPE* 

DEFINITION 

r mag : 

FT 

1 

C 

Magnitude of R^q 

r tmag 

FT 

1 

C . 

Magnitude of Ry 

r thrust > 

FT 

3 

C 

Second integral of thrust acceleration vector 
over thrusting maneuver 

r bias 

FT 

3 

C 

A position bias to account for effects of a 
rotating thrust vector 

R C1 

FT • 

3 

C 

Initial position vector input to precise 
predictor 

R C2‘ 

FT 

3 

C 

Final position output from precise predictor 

r d 

FT 

3 • 

C 

Desired thrust cutoff position 

r gd 

FT 

3 

. I 

Vehicle position vector used by guidance 

r go ' , 

FT 

3, 

C 

Position-to-be-gained including bias 
(range-to-go) 

r goxy 

FT 

3 

C 

Projection of Rqq on plane defined by Uj, and 
U Y 

r grav 

FT 

3 

C 

Second integral of gravitational acceleration 
over thrusting maneuver 
• 

*C = COMPUTED; I 

= INPUT; 0 = 
\ 

OUTPUT 





TABLE IV - DEFINITIONS OF SYMBOLS USED IN PEG-4 

(CONTINUED) 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

TYPE* 

DEFINITION 

R INT 

FT 

3 

C 

Intermediate position vector computed and used 
.in precise predictor 

Rp 

FT 

3 

C 

Predicted thrust cutoff position 

Rt 

'FT 

3 

I 

Target position in inertial coordinates 

rho mag 

NA 

I 

C 

Scalar damping factor applied to V' MISS to correct 
V G0 ' 

s 

FT 

1 

C 

Double time integral from zero to T then zero to 
. T gq of (FT/M) 

S MI$S 

NA 

1 

I 

Decimal fraction of Iv^lto be used as a con- 
vergence criterion for 1 ^mtss 1 

SA( ) 

FT 

N PHASE 

C 

Segment of S contributed by each guidance phase 

t gd 

SEC 

1 

I 

Current guidance time . 

t go 

SEC 

1 

C/O 

Time-to-go to thrust cutoff 

t goa^ J 

SEC 

N PHASE 

C 

Time-to-go until end of each phase 

t gob 

SEC 

1 

C 

Value of T gqa from previous phase 

t gop 

SEC 

1 

C 

Predicted thrust cutoff time 

*C = COMPUTED; I « 

INPUT; 0 = OUTPUT 






TABLE IV - DEFINITIONS OF SYMBOLS USED IN PEG-4 

(CONTINUED) 


FLOW CHART 


SYMBOL 

UNITS 

DIMENSION 

Tp : 

SEC 

1 

T PREV 

SEC 

1 

T RCS0FF 

SEC 

1 

V ) 

SEC 

n phase 

TT 

SEC 

1 

r 

SEC 

1 

r 

T * 

SEC 

1 

f ep . 

NA • 

3 

u T . 

NA 

3 

V 

NA 

3 


NA 

3 


NA 

3 

*C = COMPUTED; I = 

INPUT; 0 = OUTPUT 

♦ 


TYPE* 

C 

C 

I 

C 

C 

C 

C 

I 

C/O 

C 

C 

C 


DEFINITION 

Predicted time of burnout 
Previous value of Tp 

RCS cutoff time during AOA 0MS-1- parallel OMS/ 
RCS maneuver 

Estimated burn time remaining in each phase 

Time remaining until target is achieved 
Previous value of T GD 

Reference time associated with X amd X 

Unit vector in direction of Earth's axis of 
rotati on 

Resultant unit vector in the thrust direction 

Unit vector along desired thrust cutoff position 
vector 

Unit vector normal to desired trajectory plane 
(.downrange X radial) 

Unit vector along downrange (U^ X Uy) direction 
at thrust cutoff 



TABLE IV - DEFINITIONS OF SYMBOLS USED IN PEG-4 

(CONTINUED) 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

TYPE* * 

DEFINITION 

UUNIT 

NA 

3 

C 

Unit vector normal to the desired transfer 
plane in the direction of the angular 
momentum of the transfer 

V C1 

FT/SEC 

3 

C 

O 

Initial velocity input to precise predictor 

-< 

m 

X 

FT/SEC 

1 

I 

Total effective exhaust velocity 

V EXA^ ^ 

NA 

N PHASE 

C 

Effective exhaust velocity for each thrust 
phase 

V GOD 

FT/SEC 

1 

C 

Desired Vqq magnitude based upon desired cutoff 
mass for fuel depletion 

V GOMAG 

FT/SEC 

1 

C 

Magnitude of Vgg 

V GOYS 

FT 2 /SEC 2 

1 

C 

Square of out-of-plane V rn required to deplete 
OMS fuel 

V MISS 

FT/SEC 

3 

C 

Difference between predicted and desired thrust 
cutoff velocity 

7 P 

FT/SEC 

3 

c 

Predicted thrust cutoff velocity 

7 C2 

FT/SEC 

3 

c 

Final velocity output by precise predictor 

^gd 

FT/SEC 

3 

I 

Vehicle velocity vector used by guidance 

*C = COMPUTED; 

I = INPUT; 0 = OUTPUT 






TABLE IV - DEFINITIONS OF SYMBOLS USED IN PEG-4 

(CONTINUED) 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

TYPE* 

DEFINITION 

ho 

FT/SEC 

3 

C/O 

Velocity-to-be-gained, vector 

hi? 

FT/SEC 

3 

C 

Projection of Vqq vector into plane normal to 

V GRAV 

FT/SEC 

3 

C 

First integral of gravitational acceleration 
over thrusting maneuver 

Tint 

FT/SEC 

3 

C 

Intermediate velocity used by precise predictor 

V THRUST 

FT/SEC 

3 

C 

First integral of thrust acceleration vector 
over thrusting maneuver 

XLD XZ ' 

NA 

1 

I 

Fractional multiplier of circular orbital, rate 
used in computation of X xz 

A COAST 

SEC 

1 

C 

Difference between' assumed coast time and 
actual coast time 

AC-j 

FT/SEC 

1 

I 

Offset applied to true target velocity con- 
straint C-j to obtain a conic constraint 

ARj 

FT 

3 

I 

Bias applied to true target position Rj to 
achieve a conic target 

at avg 

SEC 

' 1 

C 

Integration time step for gravity prediction 

at coast 

SEC 

1 

C ■ 

Time of coast from cutoff to target 

*C = COMPUTED; 

I = INPUT; 0 = OUTPUT 






TABLE IV - DEFINITIONS OF SYMBOLS USED IN PEG-4 

(CONTINUED) 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

type* 

DEFINITION 

AT COAST(REF) . 

SEC 

1 

I 

Assumed time of coast segment 

at max 

SEC 

1 

I 

Maximum integration time step for gravity 
prediction 

at min 

SEC 

1 

I 

Minimum integration time step for gravity 
prediction 

AV S 

FT/SEC 

3 

I 

Change in accumulated sensed velocity 

e STEER . 

SEC 

1 

• I 

Convergence criterion for inhibiting steering 
output 

X . ; 

NA 

3 

C 

Unit vector in the direction of Vgg 

X 

RAD/SEC 

3 

C 

Derivative of unit vector instantaneously aligned 
along X but rotating with the desired thrust 
vector turning rate 

X MAG . 

RAD/SEC 

1 

C 

Magnitude of X • 

hi 

RAD/SEC 

1 

c 

Component of X in the X-Z plane (primer rate) 

*Y * • 

RAD/SEC 

1 

c 

Component of X along Uy 

U EARTH 

s ? 
FT /SEC 

1 

I 

Earth gravitation constant 

$max 

RAD 

I 

I 

Small angle approximation constraint 

t ( ) 

SEC 

n phase 

r 

Ratio of mass to mass flow rate for each phase 


*C = COMPUTED; I « INPUT; 0 = OUTPUT 




TABLE V 

- DEFINITIONS 

OF SYMBOLS USED 

IN PRECISE PREDICT 

FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

\ 

TYPE* 

DEFINITION 

A i. O ' 

NA 

1=1, 9 
0=1,2 

C 

A two column array used for temporary storage 
of the Legendre polynomials and the derived 
Legendre functions (latitude dependent terms) 
used in the gravity acceleration model 

ATM 

NA 

1 

C 

Flag indicating vehicle attitude is to be 
taken from a prestored table for use in drag 
and vent acceleration calculations; ATM 
always set equal to 1 

AUX 

NA 

1 

C 

Intermediate scalar variable used in gravity 
acceleration model 

r 

w i 

NA 

L = 1, 35 , 

I 

Constants which account for the tessera! 
harmonics in gravity acceleration model 

o • 

FT/SEC 2 

3 

C 

Drag acceleration vector 

DM 

NA 

1 

I 

Drag 

= 0; do not compute drag acceleration 
= 1; compute drag acceleration 

DMN 

NA 

1 

• 

C 

Intermediate variable used to compute effect 
of tesseral harmonics in gravity acceleration 
model 

[EF TO M50 (T)3 

NA 

3,3 

I 

Transformation matrix from earth-fixed .to M50 
coordinates * 

[FIFTY] 

NA 

3,3 

c 

Transformation matrix set equal to [EF TO M50 (T)] 
matrix 

F-, >F 2 ,F 3 5F 4 

NA 

1 

c 

Auxiliary variables which store the result of 
multiplying the Legendre polynomials by the 
zonal harmonic coefficients 

*C = COMPUTED; I = 

INPUT; 0 = OUTPUT 






* 


TABLE V 
♦ 

- DEFINITIONS OF SYMBOLS USED 
(CONTINUED) 

IN PRECISE PREDICT 

FLOW CHART 
SYMBOL 


UNITS 

DIMENSION 

type* 

DEFINITION 

G 

• 

FT/SEC 2 

3 

C 

Gravity acceleration vector due to the Earth's 
non-spherical shape 

^CENTRAL 


FT/SEC 2 

3 

C 

Acceleration vector due to the Earth's 
gravitational attraction as a point mass 

*F 


FT/SEC 2 • 

3 

C 

Gravitational acceleration vector at final 
state 

GMD 


NA 

1 

• I 

Flag indicating degree of gravity potential 
model desired 

= 0; spherical earth gravity acceleration 
> 0; non-spherical earth Gravity acceleration 
0 < GMD < 8 

GMO 


NA 

1 

I 

Flag indicating order of the gravity potential 
model 

0 < GMO < 8; GMO < GMD 

^PREV 


FT/SEC 2 

3 

C 

Gravitational acceleration vector at previous 
state 

,n step 


NA 

1 

C 

Number of integration steps to be used in 
PRECISE PREDICT 

R 

• 

FT 

3 

C 

Current radius vector used in ACCEL function 

*1 

r egrav 


FT 

' 3 

I 

Initial radius vector inppt to PRECISE PREDICT 


FT 

1 

I 

Earth radius magnitude used in gravity 
calculations 


*C = COMPUTED; I = INPUT; 0 = OUTPUT 



TABLE V - DEFINITIONS' OF SYMBOLS USED 

(CONTINUED) 


■ FLOW CHART 
SYMBOL 


UNITS 

DIMENSION 

TYPE* 

R F 

. 

FT 

3 

C/O 

r ef 


FT 

3 

C 

r inv 


FT 

1 

c 

RCS 


FT/SEC 2 

3 

c 

ro zero 


NA 

1 

c 

r° n 


FT/ SEC 2 

1 • 

c 


S L 

NA 

L='i , 35 

I 

T* 

SEC 

1 

C 

V 

SEC 

1 

I 

T I 

SEC 

1 

I 

U R 

NA ■ 

3 

c 

¥' 

FT/SEC 

3 

C/I 

VENT 

FT/SEC 2 

3 

c 


IN PRECISE PREDICT 


DEFINITION 

Final radius vector output by PRECISE PREDICT 

Current radius vector in Earth-fixed coordinates 

Inverse of current radius magnitude in inertial 
coordinates 

Acceleration vector due to RCS venting 

Variable used as starting value for recursive 
relations used in Pine's formulation of the 
spherical harmonics development 

Distance related term used as a starting value 
for the recursive relations used in Pine's 
formulation of the spherical harmonics 
development 

Array of constants which account for the sectorial 
harmonics in the gravity potential model 

Current time used in ACCEL function 

Final time input to PRECISE PREDICT 

Initial time input to PRECISE PREDICT 

Unit vector in direction of current radius in 
Earth fixed coordinates 

Final velocity output from PRECISE PREDICT 
and used in ACCEL 

Acceleration vector due to venting forces 





TABLE 'V 

- DEFINITIONS 

OF SYMBOLS USED 
(CONTINUED) 

IN PRECISE PREDICT 

FLOW CHART 
SYMBOL 


UNITS 

dimension 

• TYPE* 

DEFINITION 

V F 


FT/ SEC 

3 

C/0 

Final velocity output from PRECISE PREDICT 

7 I 


FT/ SEC 

3 

c/I 

Initial velocity from PRECISE PREDICT used 
in ACCEL 

VM 


nA 

1 

/ 

I 

Vent model flag 

= 0; do not compute vent acceleration 
= 1; compute vent acceleration 

zonal n 


NA 

N=1 , 8 

I 

Array of constants which account for the zonal 
harmonics in the gravity potential model 

Z IMAG X 


NA 

1=1, 9 

c 

Array which stores the effects of tessera! 
harmonics in gravity potential model. 

Z REAlj 


NA 

1=1, 9 

c 

Array which stores the effects of tessera! 
harmonics in gravity potential model 

at max 


SEC 

1 

I 

Maximum integration st'ep size 

at step 


SEC 

1 

c 

Integration step size 

^ EARTH 


ft 3 /sec 2 

1 * 

I 

Earth gravitational constant 


*C = COMPUTED; I = INPUT; 0 = OUTPUT 



TABLE VI - DEFINITIONS OF SYMBOLS USED IN LTVCON 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

TYPE* 

DEFINITION 

A 

NA 

1 

C 

Quadratic coefficient in the quadratic 

-equation " 

* 

B 

NA 

1 

C 

Linear coefficient in the V^-j quadratic equation 

b d 

NA 

1 

C 

Denominator of the coefficient function used in 
the continued fraction evaluation of the 
hypergeometric function F 

b n 

NA 

1 

C 

Numerator of the coefficient function used in 
the continued fraction evaluation of the 
hypergeometric function F 

c 

NA 

1 

C 

Constant coefficient in the V H1 quadratic equation 

C 1 ' 

FT/SEC 

f 1 

• I 

Target velocity constraint coefficient in the 
equation V RADIAL = C, + C-, V HQRIZ 

C 2 

NA 

1 

I 

Target velocity constraint coefficient in the 
equation V RABIAL =. C, + C-, V H0RIZ 

CR 

FT 

1 

C 

Magnitude of the chord connecting Rq and 

jce 

NA 

1 

C 

Cosine of the transfer angle 

D 

NA 

1 

C 

Radical of the V H1 -quadratic equation 

E 

NA 

1 • ' 

C 

Scratch variable used in the transfer time 
interval computation’ 

” F ' 

NA 

1 

c 

Hypergeometric function which is the transcendental 
portion of the normalized transfer time interval 
computation and' is evaluated by means of a 
continued fraction 


*C = COMPUTED;, I = INPUT; 0 = OUTPUT 



TABLE VI - DEFINITIONS OF SYMBOLS USED IN LTVCON 

(CONTINUED) 


FLOW CHART 

• 



. 

SYMBOL 

UNITS 

DIMENSION 

TYPE* 

DEFINITION 

f old 

NA 

1 

■ C 

Value of F at the previous level of the continued 
fraction 

I FLAG l 

’ NA 

1 

C/0 

Error flag: 0 = no error; 1 = no physical 

solution; 2 = hyperbolic trajectory 

LEVEL 

NA 

1 

C 

A measure of the level of the continued fraction 

QM 

NA 

1 

C 

Symbol representing the function subroutine 
which executes the continued fraction evaluation 
of the hypergeometric function F 

r circ 

FT 

1 

C 

Half of the semi -perimeter of the transfer 
triangle, and used as the basis for the 
normalization of the equations 

R o 

FT 

3 

I 

Initial inertial position vector 

^1 

FT 

3 

! 

Target radius vector 

R °MAG 

ft' 

1 

C 

Magnitude of Rg 

R 1 • 
'mag 

FT 

1 

c 

Magnitude of IR-j 

S . ■ * ~ 

FT ” 

1 

c 

Semi -perimeter of the transfer triangle 

S9 

■ NA 

1 

c 

Sine of the transfer angle 

U 

NA . 

1 

c 

Lambert routine independent variable 

°IN 

NA 

3 

I 

Unit vector normal to the transfer plane and - 
in the direction of the angular momentum of 
the transfer 

*N 

NA . 

3 

c 

Internal vector used to represent the unit 
normal defining the transfer plane 


*C = COMPUTED; I = INPUT; 0 = OUTPUT 



TABLE VI - DEFINITIONS OF SYMBOLS USED IN LTVCON 

(CONTINUED) 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

■ TYPE* 

DEFINITION 

%> 

NA 

* 

3 

C 

Unit vector in the direction of the inertial 
position vector R Q 

% 

NA 

3 

C . 

Unit vector in the direction of the target 
position vector R-j 

V 

NA 

1 

C 

Scratch variable used in the continued fraction 
evaluation of the hypergeometric function F. 

*0 

FT/SEC 

3 

C/0 

Required inertial velocity vector at the 
inertial position vector Rq 

Y 1 

FT/SEC 

3 

C/0 

Resulting inertial velocity vector at the 
target position vector R-j 

V CIRC 

FT/SEC 

1 

C 

Velocity magnitude of a circular orbit of radius 
R CIRC anci use< ^ ln norma ^ ization process 

V H0 

FT/SEC 

1 

C 

Normalized horizontal component of the resulting 
velocity vector V g 

V H1 

FT/SEC 

1 

C 

Normalized horizontal component of the resulting 
velocity vector V^ 

V R0 

FT/SEC 

1 

■ c- 

Normalized radial component of the required 
velocity vector Vg 

V R1 

FT/SEC 

1 

C 

Normalized radial component of the resulting 
velocity vector V-j 

W 

NA 

1 

C 

Scratch variable used in the normalized transfer 
time interval computations 

ww 

NA 

\ 

1- 

C 

Argument of the continued fraction evaluation 
of the hypergeometric function F 

*C = COMPUTED; 

I = INPUT; 0 = OUTPUT 






TABLE VI - DEFINITIONS OF SYMBOLS USED IN LTVCON 

(CONTINUED) 


FLOW CHART 
SYMBOL 

UNITS 

DIMENSION 

TYPE* 

DEFINITION 

V 

NA 

1 

C 

Scratch variable used in the continued fraction 

A * 


evaluation of the hypargeometric function F 

AT 

SEC 

1 

C/O 

Resulting transfer time interval 

A 

NA 

1 

C 

Constant parameter of the problem 

U EARTH 

FT 3 /SEC 2 

1 

I 

Earth gravitational constant 

y 

ft 3 /sec 2 

1 

• C 

Interval variable used to represent y^y^ 


*C = COMPUTED; I = INPUT; 0 = OUTPUT 
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FIGURE 2 - PEG-7 LOGIC FLOVi 


ORIGINAL PAGE IS 
OF POOR, QUALITY 
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PREDICT POSITION AND VELOCITY 
VECTORS AT CUTOFF. 


rr UPDATE POSITlOi'i AHO VELOCITY 
— VECTORS AT CUTOFF TO ACHIEVE- 
DESIRED COAST SOLUTION (8URH/ 
COAST SEQUEHCE). 



FICURE 3 - 


PEG-4 LOGIC'FLOU (PAGE 3) 






FIGURE 4 IN1T LOGIC FLOW 
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FIGURE 5 - UPVGO LOGIC FLOU 


OBIGMM} PAGE IS 
01 POOR QUALITY 
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SAVE PREVIOUS VALVE OF v 
I TIKE TO GO. 


” TEST FOR MULTI-PHASE LOGIC. 


DETERMINE SIliGLE PHASE T GQ . 




COIPUTE t AHD'SET EXHAUST 
VELOCITY FOR THIS PHASE. 



' ' W ATR 

V EXA {K PHASE' = V EX 


l - 


t = V £X /ATR 


T g0 =t 0.0 -e 




G011AG' EX\ 



TEST FOR FIRST OR SECOND PHASE 
OF GUIDANCE. 


SET FIRST PHASE BUSH TIME TO 0. 


COMPUTE FIRST PHASE BURN TIME 
AND SECOND PHASE EXHAUST VELOCITY. 


.PHASE 

< 

I NO 


LA(H) - V G0 „ AG 

T B (1) = 0.0 


V 1 ) t p.csoff ' t gd 

LA(I) = -V EX * LH (1.0 - T B 0)/tO)} 
LA(2) - V G0 ., AG - LAO) 

v exa' 2) = F QI1S / ‘W 
t( 2) = 00) - TgO)) * M/i’i 0MS 


6 ' 


FIGURE 6 - TGO LOGIC FLOtf 


0HIOTM3 vmiij is 
- OF FOOU QUALITY 















FIGURE 7 - INTEG LOGIC FLOW 
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COMPUTE CONSTANT THRUST 
INTEGRALS. 


\ 


SUM INTEGRALS OVER ALL 
.GUIDANCE PHASES. 


HAVE VALUES BEEN COMPUTED 
FOR ALL GUIDANCE PHASES? 


SA(1 } = -LA(I) * (t(I) - Tg(I)) + •'■' EXA (0 * 

OA (I) ‘ - SA(I) * LA(I) * T gqa (I)' 

QA{I) = SA( I ) * MD + W ' - 5 V EXA {1) * V 0 ’ 

* ’ t 

PA(I) = QA{I) * (t{I) * Tggg) - - 5 V exa^^ * V 1 ^ * (T GOB + 


d> 


Y — 

-HA(I) = JA(I) '* T goa O) - QAO) 

- SA(I ) ='SA(I) + L 

* t 8 (i) • 

QA(I) = QA{1) + 0 

: V 11 

* PA(I) = PA(I ) + H 

l V'i 

L = L * LA(I) 

- 

S = S + SA(I) 

- 

H = K + HA( I ) 


0*0 + 0A(1) 

• 

Q 3 Q + QA(I) 

P * P + PA(I J 





NO 



FIGURE 7 - INTEG LOGIC FLOW (PAGE Z) 


PAGE is ' 

;■ TOOK QG itn^. 


CD/3) 
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FIGUPE 7 - 1NTEG .LOGIC FLOW (PAGE 3) 


'TGIMAE PAGE IE 
- "OOR QUALITY • 
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FIGURE 8 * fURNR LOGIC FLOW 
















PREDICTOR 



FIGURE 9 - PREDICTOR LOGIC FLOW 



/• 


70 






DETERMINE PROPAGATED STATE WITH J2 
EFFECTS INCLUDED.' . 




FIGURE S - PREDICTOR LOGIC FLOW (PAGE Z) 




t 
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•' TEST FOR ABORT OR 
.1 DEORBIT MANEUVER 


ADJUST DESIRED TERMINAL 
POSITION AND TARGET VECTOR 
BIAS TO BE IN FLIGHT PLANE. • 


TEST FOR AOA 0:i$-2 OR DEORBIT — ' 
MANEUVER | 


ADJUST DESIRED TERMINAL 


. RADIUS TO ENSURE REQUIRED 
ALTITUDE 



FIGURE 10 --CORRECTOR LOGIC FLOW 
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ORIGINAL PAGE IS- 
OF POOR QUALITY 













• FIGURE 10 - CORRECTOR LOGIC FLOW (Pago 3) 
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determine: number ano size of 

INTEGRATION STEPS AND INITIALIZE 
FINAL STATE TO INITIAL STATE. 



INPUT: GKO, GKO, OK , VM , AT, 

f pp.ecisc 

R l*V T I' T f 

V PREDICT. J- 


V 

OUTPUT: Rp, Vp 

CEILING (|t; - T, 1 

/AT„, V ) © 


COMPUTE GRAVITY ACCELERATION 
FOR THE STATE AND TIME 


ACCELtGt'.D.GMO.DM , VII , ATM. ° T , Vp 


r USE "SUPER-G" ALGORITHM TO 
INTEGRATE TO THE FINAL STATE, 


T * T + AT C 


■ = R F + iT $TEP*^ V F + -5 iT STEP * G PREY* 
p = ACCEL (GMD.GMO, DM , VM . Rp, V p . T) 

F = V F + .5 * (GpREY + S F } * 6T STEP 
R = R + AT„ * * (G t - G D „rv )/ 6 



I * CEILING (A/3) IS A FUNCTION WHICH TAKES THE MAXIMUM INTEGER 
L ^ALVE_OTj A_/Uj_ ANO JNCREM_ENTS_IHAT VALVE_BYJ0NE . 


FIGURE 11 * PRECISE PREDICT LOGIC FLOW 


OB3C5S 


















FIGURE 12 - ACCEL LOGIC FLOW 


ORIGINAL PAGE IS 
OF POOR' QUALITY 
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m 


IF VENT ACCELERATION OESIRED CALL 
YEUT- MODEL 


IF DRAG ACCELERATION DESIRED CALL 
DENSITY AND DRAG F3DELS 


SU!t THE ACCELERATION VECTORS OF . 
SPHERICAL EARTH, NON -SPHERICAL * . 

EARTH, DRAG AND VENT FOP.CES 

— | Acca r g cektral + G + D ” VEHT 

‘ 

RETUR-l" - ^) 


FIGURE 12 - ACCEL LOGIC FLOW {PAGE 2) 














FIGURE 13 - ACCEL GRAY LOGIC FLOW (PAGE 2} 


ORIGINAL 

0 r p'OOB, QUALm 
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FIGURE 13 - ACCEL GRAY LOGIC FLOW (PAGE 3) 








•FIGURE 13 - ACCEL GRAY LOGIC FLOW (PAGE 4} 
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1 DETERMINE TRANSFER PLANE AND 
I SET CONSTANTS 


~i DETERMINE INITIAL AND FINAL 
j VELOCITY VECTORS 


DETERMINE TRANSFER TIME 



FIGURE 14 - "LTYCON LOGIC FLOW 


ORIGINAL PAGEJS 

• AT? POOR QUALITY 
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I 


INITIAL 


J 


DETERMINE TRANSFER PLANE, TRANSFER 
ANGLE AID SET CONSTANTS FOR 
NORMALIZING 


I FLAG = 0 

V ' l '» l 
V '!"' 1 

U RO = V R 0, 


HAG 


HAG 


• > * V R i 

C0 = B RO’.°R1 

SO = {U R0 X U R1 ) \Uj H 

CR = - R 0 { 

s = + K 


R, + CR)/2 
. U HAG ‘.HAG 


X = SIGH <SS) CR/? 


P 'CIRC " S/2 
V CIRC = v/R CIRC 


V 

^ RETURN^ 


FIGURE 15 - INITIAL LOGIC FLOW 


original PA ff Tq S 

OF POOR QUALITY 




FIGURE 16 - REqVEl LOGIC FLOW 


■SHNAE H®® 

,? POOR QUALITY 
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FIGURE IB - FUNCTION QH LOGIC FlOtf 


-.■'■aaSAE vhC&jS 


86 







PEGSUP 






LYNDO^B.«JOHN^ft*SRACE CENTER MEMORANDUM 

■ — -W- ■ — r : 


INITIATOR; 


LC-5-363 


LC/ P w "rri field: rl : 5-1 9-76 : 2091 


JM6/E. Robinson 


from; LA3/Manager for Management Integration 


Document Release Authorization 



Reqinald M. Machell 


On September 23, 1975, by memorandum LC-5-260 from this office, you were 
requested to process document release authorizations for documents generated 
by the Space Shuttle Engineering and Operations Support Contract NAS 9-13970 
according to the following criteria: 

a. "May be announced in STAR (or CSTAR if classified)." 

b. "May be made publicly available, except (1) classified material which 
should be released to security qualified requestors only, and (2) documents 
which may be released overseas that come under the jurisdiction of the Federal 
Early Domestic Dissemination Program." 

The above referenced contract has been superseded by NAS 9-14960, but the 
processing of document release authorizations should continue to be guided by 
the same criteria. If any additional clarification of this criteria is needed, 
such requests should be addressed to LC/R. B. Merrifield, extension 2091. 


JSC Form 3180 (Mar 75) 


INCREASED PRODUCTIVITY ~ LOWER COST 


PAGE 1 Of 










